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ABSTRACT 

We evaluate the minimum energy input rate that starbursts require for expelling their newly 
processed matter from their host galaxies. Special attention is given to the pressure caused 
by the environment in which a galaxy is situated, as well as to the intrinsic rotation of the 
gaseous component. We account for these factors and for a massive dark matter distribution, 
and develop a self-consistent solution for the interstellar matter gas distribution. Our results are 
in excellent agreement with the results of Mac Low & Ferrara (1999) for galaxies with a flattened 
disk- like ISM density distribution and a low intergalactic gas pressure (PiGAt/k < 1 cm -3 K). 
However, our solution also requires a much larger energy input rate threshold when one takes 
into consideration both a larger intergalactic pressure and the possible existence of a low-density, 
non-rotating, extended gaseous halo component. 

Subject headings: ISM: bubbles, ISM: abundances, general - starburst galaxies 



1. Introduction 

The observational evidence of powerful star- 
bursts in dwarf galaxies has led to the idea that, 
due to their rather shallow gravitational poten- 
tial, supernova (SN) products and, even the whole 
of the interstellar medium (ISM), may be easily 
ejected from the host dwarf systems, causing the 
contamination of the intra-cluster medium. This 
issue has been addressed in several papers. Do 
Young & Gallagher (1990) and De Young & Heck- 
man (1994) concluded that typically 1/6 of the 
shock-wave's swept-up mass may escape its host 
galaxy, and coherent SN explosions may even to- 
tally destroy a galaxy's ISM. More recently, how- 
ever, Silich & Tenorio-Tagle (1998) and D'Ercole 
& Brighenti (1999) have shown that within a dark 
matter (DM) and a low-density extended gaseous 
halo, a galactic wind does not easily develop, mak- 
ing the removal of the galaxy's ISM very ineffi- 
cient. Recently, Mac Low and Ferrara (1999, here- 
after MF) have used several observational rela- 
tions to define the galactic dark matter component 
by means of the central DM density (p c ), its maxi- 
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mum extent (Rh) and its characteristic scale (R c )- 
They made quantitative conclusions regarding the 
energetics required for the ejection of metals and 
of the ISM ("blow-away") from galaxies with an 
ISM mass in the range {M ISM = 10 6 - 10 10 M ). 

In none of the above mentioned studies, how- 
ever, have the roles of the extragalactic environ- 
ment nor of the rotation of the interstellar gas 
been examined in full detail. Meanwhile, Babul 
& Rees (1992) have indicated that the extragalac- 
tic medium pressure (PiGM/k) may vary between 
lO^ 1 and 10 4 cm~ 3 K, making it one of the key 
parameters in the galactic-extragalactic medium 
connection. 

In the present paper we reexamine the effi- 
ciency of starbursts in ejecting their newly pro- 
duced metals, and in ejecting their whole ISM, 
from dwarf and normal galaxies, taking into ac- 
count new boundary and initial conditions. These 
follow from the pressure equilibrium condition at 
the galactic-extragalactic boundary, and the rela- 
tive ability of rotation and interstellar gas pressure 
to balance different fractions of the galaxy's grav- 
ity. In particular, we show that in a self-consistent 
model, the ISM's gas velocity dispersion and its 
density at the galaxy center are no longer free pa- 
rameters, but are defined by the total ISM mass 



(Mism), the intergalactic pressure (Pigm), and 
the fraction of the radial component of gravity bal- 
anced by pressure gradients. It follows from our 
considerations that the models of MF which have 
a low velocity dispersion and a flattened ISM disk- 
like density distribution, represent a limiting case 
in the more general scenario of starburst dwarf 
galaxy evolution. 

Section 2 presents the derivation of our model 
galaxy. Section 3 is devoted to a search for the 
threshold energy input rate required to expel ei- 
ther metals or the ISM from galaxies with a total 
ISM mass in the range of 10 6 to 10 9 M . Section 
4 is a discussion of our results. 

2. The Galaxy Model 

Our initial model (Silich & Tcnorio-Tagle, 1998; 
hereafter referred to as Paper I) included stars, 
dark matter, and several isothermal interstellar 
gas components related to the ionized, neutral and 
molecular ISM phases. It was in fact an exten- 
sion to the model developed by Morita (1982), 
Tomisaka & Ikeuchi (1988), Tomisaka & Bregman 
(1993) and Suchkov et al. (1994). To be consis- 
tent with the approach of MF we have simplified 
our model by using only the gravitational field im- 
posed by the dark matter (DM) component and 
using only one of the ISM phases with an effective 
velocity dispersion Gism- 

The spheroidal DM mass distribution can be 
approximated by 



Pdm(t) 



1 + {r/Rc) 2 ' 



(1) 



The DM halo radius (Rh), characteristic scale 
height (R c ) and central density (p c ) are given by 
MF: 
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where H is the Hubble constant in 100 km 
s _1 units. We adopt H=0.65 throughout the pa- 
per. 



On the other hand, the interstellar gas density 
and pressure distributions are then defined by (see 
paper I) 
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where p is the ISM density at the galactic center, 
V h = ^J 2G r" m is the escape velocity at the DM 
boundary Rh, G is the gravitational constant and 
Mdm is the DM mass, which is given by 
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X = F(u) - e 2 F(r) - (1 - e 2 )^(0), (8) 



where r = \J (x 2 + y 2 ) and w = \J (x 2 + y 2 + z 2 ) 
are the cylindrical and spherical radii, respec- 
tively, and the dimensionless radii = uo/R c and 
yh = Rh/ R c - The factor (1 — e 2 ) is the fraction of 
the radial component of gravity that is balanced 
by pressure gradients (Tomisaka & Ikeuchi, 1988). 
The case e = 1 implies a full balance between the 
radial component of gravity and the circular rota- 
tion of the galaxy, whereas e = implies spher- 
ically symmetric systems without rotation. The 
function F was defined in paper I and, in the case 
of a potential exerted by DM only, it becomes: 
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The DM mass is related to the visible mass by 
(MF) 



M DM = 3.47 x 10* 
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One can also include a transition from a rotat- 
ing disk to a non-rotating spherical halo by ex- 
pressing the e in equation (8) as 



e = e c /e^[{z/H z f + {r/H r f 
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where e c is the value of e at galactic center. This 
prevents the funnel-like density distribution that 
develops around the symmetry axis (see Suchkov 
at al. 1994). However, to be consistent with the 
considerations of MF, we have set the scale- lengths 
H z and H r to be much greater than or comparable 
to the DM halo radius. 

Following MF we assume that the interstellar 
gas extends out to a cutoff radius Rism given by 
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However, we further assume that the galaxy is sta- 
ble and that the ISM is in pressure balance with 
the intercluster gas at the galaxy outer boundary 
(Pism = Pigm)- This defines the ISM gas density 
at the galaxy outer edge to be 
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Thus, given the total mass of the ISM (Mjsm) 
and the pressure of the intergalactic medium 
(Pigm), only po and Cism must be found for 
fully defining the ISM mass distribution. This can 
be achieved by solving a system of two nonlinear 
equations 

/>Rh rRh 
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Jo Jo 

Pism(Rism,0) = Pigm, (16) 

where pism{t, z) and Pism{t, z) are defined by 
the equations (5) and (6). Note that Pism{t,z) 
should be set equal to zero whenever its calculated 
value falls below pq. The equations can be solved 
by an iterative method of high accuracy. Thus, in 
a self-consistent model, the ISM gas velocity dis- 
persion and gas density at the galaxy's center are 
no longer free parameters, but are determined by 
the total ISM mass Mjsm, the extragalactic pres- 
sure Pigm, and the fraction of the radial compo- 
nent of gravity balanced by pressure gradients. 

Figure 1 presents as an example, the ISM den- 
sity distribution for Mjsm — 10 8 M Q models with 



different fractions (e 2 ) of the radial component of 
gravity balanced by the centrifugal forces. The fig- 
ure shows the smooth transition from a flattened 
disk-like system dominated by fast rotation (Panel 
a) to a spherical, non-rotating galaxy (Panel c). 

3. Energy requirements 

Once a nuclear coeval starburst has all its stars 
on the main sequence, it begins to build a su- 
perbubble. In a constant density medium, the 
shock enters a phase of continuous deceleration 
from the very beginning of the superbubble evolu- 
tion. However, if the starburst is placed in a plane 
with a stratified decreasing gas density, it will ex- 
pand furthest in the direction of least resistance, 
deforming the originally round superbubble into 
an elongated remnant along the symmetry axis. 
In the case of an exponentially-decreasing density 
distribution, the remnant will then break out if 
the leading shock reaches a couple of galactic scale 
heights (H e ) while maintaining a supersonic speed, 
V s (Kompaneets 1960, Koo & McKee 1992), while 
in a Gaussian atmosphere the shock acceleration 
starts once it reaches a galactic scale-height H g 
(Koo & McKee 1990). As Vs = {Pbubbie/pz) 05 , 
where PbubUe is the interior pressure, and pz is 
the ISM density that decreases sharply along the 
Z axis, the shock is forced to accelerate. Once 
in the acceleration phase, neither the external gas 
pressure nor the galaxy's gravitational potential 
can restrain the shock, which ultimately will ac- 
quire speeds larger than the galactic escape ve- 
locity (V esc ). This occurs even if V s < V esc prior 
to the acceleration phase that begins at Z ~ the 
gaseous scale-height. It is also well known that the 
swept up matter immediately behind the acceler- 
ating shock, also accelerates, producing the on- 
set of Rayleigh- Taylor instabilities that fragment 
the swept-up matter shell. The remnant is then 
able to "blowout", driving its high pressure gas 
between the fragments into the halo of the host 
galaxy, where it once again will push the shock 
(see Tenorio-Tagle & Bodenheimer 1988, and ref- 
erences therein). 

The minimum starburst energy input rate 
(L cr ), that leads to breakout, continuous accel- 
eration, and expulsion of matter from galaxies 
into the IGM, either from a Gaussian or an expo- 
nential atmosphere was studied by Koo & McKee 
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(1992). The authors used the threshold luminosity 
criterion (Lbreakout) which assures that the shock 
will reach one scale height H while maintaining a 
supersonic speed: 

Lbreakout = 17.9p H 2 C] SM ergs s" 1 (17) 

For an exponential atmosphere, Koo & McKce 
(1992) demand an L cr > iLbreakout- In a Gaussian 
density distribution the acceleration phase starts 
earlier, and thus we shall assume L cr > Lbreakout- 
Note, that in the flattened disk-like configurations, 
the characteristic scale height H may exceed the 
thickness of the galaxy's disk Zism, requiring that 
H in the equation (17) be replaced by Zism- 

3.1. The Density Distribution along the 
Axis of Symmetry 

Applying the breakout criterion to our model 
galaxies required accurate fits to the resultant den- 
sity distributions (see Section 2) along the symme- 
try axis. 

3.1.1. Disk-like density distributions 

In the case of flattened disk-like distributions, 
one can show that near the galaxy's plane the re- 
sultant gas density distribution along the symme- 
try axis is fitted well by a Gaussian function: 

Pism = p exp[-(Z/H g )% (18) 

with the characteristic scale height 



the following two equations: 
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Further away from the galaxy's plane, the density 
distribution can be well approximated by an ex- 
ponential function 



Pism = Pg exp[(Z /S M - Z)/H e 
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where Zism is the disk thickness along the Z axis. 
As the match to the density distribution requires 
two different functions, we have restricted the pa- 
rameter H e by the condition that the density gra- 
dient, as given by the functions (18) and (20), 
must be smooth and continuous at the crossing 
point (Zf). This condition can be represented by 
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Solving for H e and Zf gives: 
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which uniquely define the intersection point Zf 
and the fitting parameter H e required in estimat- 
ing the minimum energy input rate for breakout. 

The ISM distribution along the Z axis for dif- 
ferent disk-like galaxies and different extragalac- 
tic gas pressures (solid lines) and their analytic fit 
(dashed lines) are shown in the Figures 2 a and b. 

3.1.2. Spherical density distributions 

The ISM density distribution for different 
galaxies under the assumption of e = 0, and dif- 
ferent extragalactic gas pressures are shown in 
Figures 2 c and d. In these cases the density dis- 
tribution far from the galaxy's plane cannot be 
fitted by a Gaussian or an an exponential func- 
tion. Thus to estimate the minimum energy input 
rate required to eject the metals from these galax- 
ies we have only used the numerical scheme of 
paper I. In all cases we have assumed a constant 
mechanical energy input rate during the starburst 
supernova phase (4 xlO 7 yr). 

3.1.3. The energy requirements 

Our galaxy models are summarized in Table 
1. The various models are labeled with several 
indexes that represent the logarithm of the ISM 
mass considered, e values that range between 
and 0.9, and intergalactic gas pressure (PiGM/k) 
values that arc assumed to be cither 1 or 100 cm -3 
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K. In the Table, column 1 gives the model identi- 
fication. Columns 2 and 3 list the DM and ISM 
masses, column 4 gives the ISM cutoff radius, col- 
umn 5 the intra-cluster gas pressure (Pigm), and 
columns 6 and 7 show the interstellar gas veloc- 
ity dispersion and escape velocity at the galaxy 
boundary Rism- The central gas number density 
is shown in column 8 and column 9 shows the ISM 
extent along Z axis (Zjsm)- 

Figure 3 shows our energy estimates resultant 
from the numerical integration of the hydrody- 
namic equations for a coeval starburst with a con- 
stant energy input rate during the first 4 x 10 7 yr 
of the evolution, and the initial density distribu- 
tions shown in Figure 2 a - d. Figure 3 considers a 
range of galactic ISMs (from 10 6 to 10 9 M Q ) and 
extreme values of e (= 0.9 for flattened disk-like 
density distributions and = for spherical galax- 
ies without rotation) and also for a range of val- 
ues of Pigm Ik (equal to 1 and 100 cm^ 3 K). Our 
results for disk-like gas distributions are in excel- 
lent agreement with the findings of MF, particu- 
larly for the case of Pigm I k = 1 cm~ 3 K. Note, 
however, that the threshold energy input rate in 
our models depends on the intergalactic gas pres- 
sure, Pigm, and becomes rapidly larger for smaller 
Mism and larger Pigm- 

In the case of the flattened disk-like configura- 
tions (e = 0.9; see Table 1), our numerical exper- 
iments indicate that the shock acceleration phase 
starts after 1.2 - 1.4 H g , and thus it is the Gaus- 
sian rather than the exponential outer part of the 
disk-like density distributions the ones that define 
the onset of shock acceleration. Our results are 
in good agreement with the analytic criterion that 
ensures breakout (equation (17)) from the density 
distributions shown in Figure 2 a and b, particu- 
larly when the values of the Gaussian scale-height 
(H = H g ) are used (filled symbols in Figure 3). 

In Figure 3, values below each line imply total 
retention, while the region above each line indi- 
cates the expulsion of the hot superbubble inte- 
rior gas (the new metals) out of disk-like distri- 
butions (e = 0.9), and of the new metals and the 
whole of the ISM in the spherical (e = 0) cases. 
Each line that separates the two regions marks 
the minimum energy input rate needed from a co- 
eval starburst to reach the outer boundary of a 
given galaxy, regardless of the time that the rem- 
nant may require to do so. In the case of disk- 



like configurations, this energy input rate warrants 
breakout and the ensuing continuous acceleration. 
In most of these cases the ejection of the hot gas 
into the intergalactic medium occurs before the 
starburst supernova phase is over (4 xlO 7 yr). In 
galaxies with e = 0, however, the evolution time 
considered largely exceeds the starburst supernova 
phase. This causes the decelerating shell of inter- 
stellar swept up matter to reach the galaxy edge 
by means of the momentum gathered during its 
early evolution. Figure 4 shows the results for a 
Mism = 10 8 M galaxy under two extremes of in- 
tergalactic pressures [Pigm I k = 1 and 100 cm~ 3 
K) as they become expelled by the lower limit me- 
chanical energy input rate displayed in Figure 3. 
Throughout its evolution, the remnant in the high 
pressure galaxy undergoes a continuous decelera- 
tion, while the low pressure system is able to expe- 
rience the acceleration promoted by breakout. In 
both cases however, the time needed to reach the 
galaxy outer boundary is much larger than the su- 
pernova phase and thus both remnants slow down 
by conservation of their own momentum, to reach 
the galaxy boundary and experience total disrup- 
tion after 100 Myr and 200 Myr of evolution in the 
high and low intergalactic pressure cases, respec- 
tively. 

The right-hand axis in our summary Figure 3 
indicates the mass of the coeval starburst able to 
produce the energy input rate indicated in the left- 
hand axis. This has been scaled from the synthesis 
models of Leitherer & Hcckman (1995) for a 10 6 
M Q starburst under the assumption of a Salpeter 
IMF with a lower mass limit of 1 M Q , an upper 
stellar mass equal to 100 M Q , and a mctalicity Z = 
0.1. In this case, the log of the mechanical energy 
input rate in erg s _1 is 40.5, and scales linearly 
with the starburst total mass. 

In Figure 3, we have also indicated the results of 
another two calculations where we assumed a grav- 
itational potential corresponding to a Mism = 10 9 
M galaxy. In these two models however, we have 
adopted a spherical halo with total masses of 10 8 
M Q and 10 7 M Q and an intergalactic gas pres- 
sure equal to Pigm I ^ — 1 cm -3 K. The derived 
minimum energy input rate (or minimum mass of 
the coeval starburst required to expel the ISM) 
clearly indicates that it is the mass of the halo, 
instead of the disk component, what determines 
the minimum energy input rate necessary for ejec- 
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tion. Figure 3 shows that a halo with only 1/100 
of the mass of the disk-like configuration (Mism 
= 10 9 Mq) increases the value of the minimum 
energy input rate required for mass ejection by 
approximately an order of magnitude. An energy 
input rate nearly two orders of magnitude higher 
is needed if the spherical halo has a mass of 1/10 
of the disk-like configuration. 

4. Discussion 

We have described a simple, self-consistent, ap- 
proach to build the ISM gas distribution of galax- 
ies, accounting for the gravitational potential ex- 
erted by a massive DM component, as well as for 
rotation of the ISM and the intergalactic gas pres- 
sure. The numerical experiments and analytical 
estimates have then showed that the final fate of 
the matter ejected from a starburst region, as well 
as that of the shocked ISM, is highly dependent on 
the boundary conditions. Following a continuous 
transition from a fast rotating, thin, disk-like to a 
spherical non-rotating interstellar gas distribution, 
we have found that the ISM to be more resistant 
to ejection than estimated in earlier papers. 

We have shown (see Figure 3) that superbub- 
bles evolving in galaxies that have a gaseous disk- 
like density distribution are likely to undergo the 
phenomenon of breakout. This allows them to 
accelerate and expel all of their newly produced 
metals, and perhaps even a small fraction of the 
interstellar medium, into intergalactic space. On 
the other hand, much larger energy input rates, or 
more massive coeval starbursts (up to 3 orders of 
magnitude larger), are required to provoke break- 
out or push a shell to the galaxy outer bound- 
ary for a spherically-symmetric ISM mass distri- 
bution. Even low mass (~ 1% - 10% of the to- 
tal ISM mass), non-rotating subsystems increase 
the energy requirements by more than an order 
of magnitude. This makes the low density haloes, 
rather than DM itself, the key component in the 
evolution of dwarf galaxies. 

The halo properties, on the other hand, are 
highly dependent on the gravitational potential. 
The presence of an extended DM component 
makes the halo mass distribution more smooth and 
extended. However, the halo parameters are also 
dependent on the properties of the intergalactic 
medium. To withstand a high intergalactic pres- 



sure, the halo has to support itself with a high 
random motion, which leads to a more homoge- 
neous gas distribution compared to that found in 
low pressure surroundings. 

Clearly the energy input rates derived here are 
lower limits to the amounts required for expelling 
matter from a galaxy. Particularly because only 
one component of the ISM was considered and be- 
cause the central densities adopted are well be- 
low the values expected for the star forming cloud 
where the starburst originated. Our estimates 
thus neglect the effect of the starburst plowing 
into the parental cloud material. These are lower 
limits also because we adopted a constant energy 
input rate (see Strickland & Stevens 2000) and be- 
cause our approach neglects an additional cooling 
by mass-loading process (Hartquist et al. 1986) 
and the presence of a magnetic field which also 
could inhibit expansion (Tomisaka 1998). 

The indisputable presence of metals (in what- 
ever abundance) in galaxies implies that the su- 
pernova products cannot be lost in all cases. Note 
in particular that many well known disk galaxies 
have a high metal abundance and a large num- 
ber of centers of star formation. Most of these 
exciting star clusters are more massive than the 
10 4 M Q lower limit established by MF and the 
current paper as the minimum starburst mass re- 
quired to cause mass ejection in the case of disk- 
like systems. This lower limit for disk-like galax- 
ies with Mism < 10 9 M© (see Figure 3) implies 
that starbursts even smaller than the Orion cluster 
would break through the galaxy outer boundary 
and eject their supernova products into the inter- 
galactic medium. Nevertheless, disk-like galaxies 
can avoid losing all their freshly produced metals 
by having a halo component, neglected in former 
studies, that acts as the barrier to the loss of the 
new metals. 

The haloes, despite acting as the barrier to the 
loss of the new metal, have rather low densities 
(< n >~ 10~ 3 cm~ 3 ) and thus have a long re- 
combination time (t rec — 1/ '{arihaio)', where a is 
the recombination coefficient) that can easily ex- 
ceed the life time of the HII region (tnn = 10 7 
yr) produced by the starburst. In such a case, the 
haloes may remain undetected at radio and op- 
tical frequencies (see Tenorio-Tagle et al. 1999), 
until large volumes are collected into the expand- 
ing supershells. Note that the continuous shape 
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that supershclls present in a number of galaxies, 
while remaining attached to the central starburst, 
and their small expansion velocity (comparable 
or smaller than the escape velocity of their host 
galaxy) imply that the mechanical energy of the 
star cluster is plowing into a continuous as yet un- 
detected medium. Other authors (see MF 1999) 
have argued that blowout could leave fragments 
of the dense shell behind and that these observed 
with poor resolution may appear to form a con- 
tinuous shell much smaller than the true extent of 
the shocked region. Clearly, X-ray observations of 
these supershells will help to decide if the hot pro- 
cessed material is enclosed or not by the expanding 
shell. 

Thus the answer, the true limit for mass ejec- 
tion from galaxies, must lie between the two ex- 
treme cases that we have investigated here. Note 
however, that in the presence of a halo, it is the 
mass of the halo that sets the limiting energy in- 
put rate required for mass ejection, and not the 
mass of the disk-like component. This argument 
applies to all galaxies whether spirals, amorphous 
irregulars or dwarfs. 
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Fig. 1. — The model galaxy. Logarithmic contours showing the density distribution for a galaxy with a 
Mism = 10 8 M Q and a PiGM/k = 1 cm~ 3 K and several values of e (= 0.9, 0.5, and 0.0). The latter 
correspond to the fraction / =1 e 2 equal to 19%, 75% and 100% of the radial component of gravity balanced 
by pressure gradients in panels a, b and c, respectively. 
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Fig. 2. — The density distribution. Solid lines in Panels a and b show the density distribution along the 
symmetry axis, derived for disk-like galaxies (see Table 1) with e = 0.9 and a mass Mism = 10 6 M Q - 10 9 
Mq, and intergalactic pressures equal to 1 and 100 cm~ 3 K, respectively. The dashed lines represent the 
corresponding fit to the density distributions. Panels c and d show the corresponding density distributions 
that resulted for galaxies without rotation (e = 0). 
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Fig. 3. — Energy estimates. The log of the critical mechanical luminosity, and of the starburst mass, required 
to eject matter from galaxies with a Mism in the range 10 6 - 10 9 M©. The lower limit estimates are shown 
for galaxies with extreme values of e (= and 0.9) and for two values of the intergalactic pressure Pigm /k 
= 1 cm~ 3 K (solid lines) and Pigm I k — 100 cm~ 3 K (dashed lines). The resolution of our numerical search 
is AlogL cr = 0.1. Each line should be considered separately as they divide the parameter space into two 
distinct regions: a region of no mass loss that is found below the line and a region in which blowout and 
mass ejection occur that is found above the line. Also indicated on the right-hand axis are the energy input 
rates required for a remnant to reach the outer boundary of a halo with mass 10 8 M© and one with mass 
10 7 M©, for the case of a gravitational potential provided by Mdm = 9.1 xlO 9 M© that can hold an Mism 
= 10 9 M©. The filled squares and triangles represent the analytical energy input rate estimated by means of 
equation (17), using the H g values derived for the central Gaussian part of the disk-like density distributions 
shown in Figure 2 a, b. 
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Fig. 4. — The maximum expansion speed measured along the symmetry axis Z for 10 8 M Q galaxy with 
spherical ISM density distribution. The solid line indicates the expansion velocity for P/GM/k=100, whereas 
the dotted line represents the shell expansion for the P/GAf/k=l model. The arrow indicates the end of the 
SB supernova phase. Note that the remnants require a much longer time to reach the outskirts of their host 
galaxy. 
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Table 1: Galaxy parameters. 
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